Hybrid protoneutron stars with the Dyson-Schwinger quark model 
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We study the hadron-quark phase transition at finite temperature in the interior of protoneutron stars, com- 
bining the Dyson-Schwinger model for quark matter with the Brueckner-Hartree-Fock approach for hadronic 
matter. We discuss the dependence of the results on different nuclear three-body forces and on details of the 
quark model. A maximum mass exceeding two solar masses can be obtained with a strong three-body force 
and suitable parameter values in the Dyson-Schwinger model. With a hybrid configuration, the maximum mass 
of protoneutron stars is larger then that of cold neutron stars, such that a delayed collapse might be possible in 
principle. 
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I. INTRODUCTION 

It is generally believed that a neutron star (NS) is formed 
as a result of the gravitational collapse of a massive star 
(M > 8M0) in a type-II supernova |T, v||. Just after the core 
bounce, a protoneutron star (PNS) is formed, a very hot and 
lepton-rich object, where neutrinos are temporarily trapped. 
The following evolution of the PNS is dominated by neu- 
trino diffusion, which results first in deleptonization and sub- 
sequently in cooling. The star stabilizes at practically zero 
temperature, and no trapped neutrinos are left 031 H . 

The dynamical transformation of a PNS into a NS could be 
strongly influenced by a hadron-quark phase transition in the 
central region of the star [5,-^,]. Calculations of PNS structure, 
based on a microscopic nucleonic equation of state (EOS), 
indicate that for the heaviest PNS, close to the maximum mass 
(about two solar masses), the central baryon density reaches 
values larger than 1/fm . In this density range the nucleon 
cores (dimension w 0.5 fm) start to touch each other, and it is 
likely that quark degrees of freedom will play a role. 

In previous articles ll9l- [llll we have studied static prop- 
erties of PNS using a finite-temperature hadronic EOS in- 
cluding also hyperons iJ^ derived within the Brueckner- 
Bethe-Goldstone (BBG) theory of nuclear matter. An even- 
tual hadron-quark phase transition was modeled within an ex- 
tended MIT bag model 0131 [Till . In the present work we con- 
sider a more sophisticated quark model, the Dyson-Schwinger 
model (DSM)tI5.-Jiy. 

The study of hybrid stars is also important from another 
point of view: Purely nucleonic EOS are able to accommodate 
fairly large (P)NS maximum masses l.20.-,24.1 . but the appear- 
ance of hyperons in beta-stable matter could strongly reduce 
this value [24-28]. In this case the presence of non-baryonic, 
i.e., "quark" matter would be a possible manner to stiffen the 
EOS and reach large NS masses. Heavy NS thus would be 
hybrid quark stars. 

The paper is organized as follows. In section|II]we review 
the determination of the baryonic EOS in the BHF approach 
at finite temperature. Section [III] concerns the QM EOS ac- 
cording to the DSM, comparing also with the MIT bag model 
for reference. In section |IV] we present the results regarding 
(P)NS structure, combining the baryonic and QM EOS for 



beta-stable nuclear matter via a Gibbs phase transition con- 
struction. SectionlVlcontains our conclusions. 



II. EOS OF NUCLEAR MATTER AT FINITE 
TEMPERATURE IN THE BRUECKNER APPROACH 

The formalism proposed by Bloch & De Dominicis [29] in 
order to solve the nuclear many-body problem at finite tem- 
perature is the closest one to the Brueckner-Bethe-Goldstone 
(BBG) expansion [3U]. In this approach the essential ingredi- 
ent is the two-body scattering matrix K, which along with the 
single-particle potential U satisfies the self-consistent equa- 
tions 
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where / = 1,2,... generally denote momentum, spin, and 
hadron species. Here V is the two-body interaction, W ^Ei + 
E2 represents the starting energy, £,■ — kf/2mi + Ui{ki) the 
single-particle energy, and n{i) is a Fermi distribution at finite 
temperature. For given partial densities of all hadron species 
pi {i = «,p, A,E, . . .) and temperature T = l/j3, Eqs. ([T]) and 
(|2|l have to be solved self-consistently along with the follow- 
ing equations for the auxiliary chemical potentials /!,, 
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The grand-canonical potential density (O in the Bloch-De 
Dominicis framework can be written as the sum of a mean- 
field term and a correlation contribution [30|, 

C0 = ~L{L ^in[l+e-P^^--f'-^)+n,{k)U,{k) (4) 

+ y ^./^(2A,-w) Tr^ (^arctan [jf{W)n5{Ho - IV)] ) | , 



where the trace Tr2 is taken in the space of antisymmetrized 
two-body states and the two-body scattering matrix J€ is de- 
fined by 
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By keeping only the first term in the expansion of the arctan 
function |30], and neglecting a series of terms proportional to 
n(fe) [1 ~n(k)\ (or powers of it), which turn out to be negligi- 
ble in the temperature and density ranges relevant for (P)NSs, 
then the correlation term simplifies and reduces to 
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which are fitted to the available experimental NY scattering 
data 11341 [35i1 . Those potentials have been widely used in the 
BBG zero-temperature calculations 11241 l25l Ba BTIl . It turns 
out that at zero temperature only A and Yr hyperons appear 
in the NS matter up to very large densities. We therefore also 
restrict the present study to these two hyperon species, ne- 
glecting the appearance of thermal Y? and E+ and heavier hy- 
perons. We remind the reader that hyperonic two- and three- 
body forces are neglected in our calculations, due to the lack 
of experimental data. 

In practice we employ empirical analytical parametriza- 
tions of the free energy density, Eq. (|2), as a function of partial 
densities and temperature, as detailed in Refs. lUU [l2l [l4ll . 



which defines the grand-canonical potential in total analogy 
with the BBG binding potential, just using Fermi functions 
instead of the usual step functions at zero temperature. 
In this framework, the free energy density is then 
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and all other thermodynamic quantities of interest can be com- 
puted from it, namely the "true" chemical potential /x,-, pres- 
sure p, entropy density s, and internal energy density e read 
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To solve Eq. ([T]i one needs as input the interactions between 
hadrons. In this paper we have used as nucleon-nucleon two- 
body force the Argonne Vig potential 131]. As widely dis- 
cussed in the literature, in order to reproduce correctly the 
nuclear matter saturation point po ~ 0.17 fm^^, E/A « — 16 
MeV, three-body forces (TBF) among nucleons are usually 
introduced. In the BBG approach the TBF are reduced to 
density-dependent two-body forces by averaging over the po- 
sition of the third particle, assuming that the probability of 
having two particles at a given distance is reduced accord- 
ing to the two-body correlation function. One should notice 
that both two- and three-body forces should be consistent, i.e., 
use the same microscopic parameters in their construction. 
In this paper, we will discuss results obtained following this 
approach to TBF ll32ll . which are based on the same meson- 
exchange parameters as the two-body potential. For complete- 
ness, we also present results based on phenomenological mod- 
els for TBF, which are often adopted. In particular, we use the 
phenomenological Urbana TBF |33], here denoted by UTBF 
(UIX), which consists of an attractive term due to two-pion 
exchange with excitation of an intermediate A resonance, and 
a repulsive phenomenological central term. 

As far as the hyperonic sector is concerned, we employ the 
Nijmegen soft-core nucleon-hyperon (NY) potentials NSC89, 



III. QUARK MATTER 
A. Dyson-Schwinger equations approach 

For the deconfined quark phase, we adopt a model based 
on the Dyson-Schwinger equations of QCD, which provides a 
continuum approach to QCD that can simultaneously address 
both confinement and dynamical chiral symmetry breaking 
Ill5l[l6ll . It has been applied with success to hadron physics in 
vacuum Ill7ll38l - l40ll . and to QCD at nonzero chemical poten- 
tial and temperature [18, 4ll-l44ll. Recently efforts have been 
made to calculate the EOS for cold quark matter and compact 
stars 1 19, 45, 46], and in this paper we extend the work of 
Ref. 1 19] to the finite-temperature case. 

Our starting point is QCD's gap equation for the quark 
propagator S{p, co„;ii,T) at finite quark chemical potential jj. 
and temperature T, which in the imaginary-time formalism 
reads ' 

S{p,co„;n,Ty^ = np + iyA{oi„ + ipi) +m^ + l.{p,Oin\pi,T) 

(12) 
with the self-energy expressed as 



Y{p,COn-llJ) = T Y, 



d q 

,=_^ (2^ 



g^Dpaik;ii,T) (13) 



x—YpS{q;^J)K{q,p;^,T), 

where ©„ ~ (2n + l)nT is the Matsubara fre- 
quency for quarks, p = {p,co„), q = {q,C0i), and 
k — {k, (0„_i) — {p — q,(0„ — (Oi). Here, g is the cou- 
pling strength, Dpa{k;iJ.,T) is the dressed gluon propagator, 
and r'^{q,p;iJ.,T) is the dressed quark-gluon vertex. 
Moreover, X" are the Gell-Mann matrices, and m^ is the 
current-quark mass. As we employ an ultraviolet-finite 
model, renormalization is unnecessary. 

The kernel, Eq. ( fTlT l, depends on the gluon propagator and 
the quark-gluon vertex. Little is known about them at T ^Q 



In our Euclidean metric: {Yp,Ya} = 2Spa; Yp = 7p' Ys = 74717273' ah = 
Y,i=[ Oibc ab = Y.i=[ Cibi', and Pp timelike =^ P^ < 0. 



and jU 7^ 0, hence we follow the DSM that has been success- 
fully applied to hadron physics at jJ. = T = 0. For the quark- 
gluon vertex we use the "rainbow approximation" 



r''„{q,p-^,T)^—ya., 
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which is the first term in a symmetry-preserving scheme to 
calculate meson properties [47]. For the gluon propagator 
in vacuum, we consider a Gaussian-type effective interaction 
ll48[l in Landau gauge 
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with the transverse projector Ppa{k) = 5pa —kpka/k^- It 
is a finite-width representation of the Munczek-Nemirovsky 
model [49] used in Ref. [46], and expresses the long-range 
behavior of the renormalization-group-improved effective in- 
teraction in Refs. 1.17. .50. .51.] . Generally, at finite chemical 
potentials and temperatures, the gluon propagator can be de- 
composed into electric and magnetic parts, 

g^D^y{k,a„)^Pj,^DTik\D.l)+Pl;^DLik\D.l) (16) 

with the gluon Matsubara frequency i2„ ~2nnT, 
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P^y. At jU = and low temperatures, there 
are indications [52] that the gluon propagator is insensitive 
to temperature. In the interior of PNSs, the typical chemical 
potential is a few hundred MeV, which is much larger than the 
typical temperature, a few tens MeV. Therefore, we use the 
extended form |il_9|] of Eq. (flST l 
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with a quantifying the asymptotic freedom at large chem- 
ical potential but without any temperature effects. As in 
Ref. 11191] . we choose the set of parameters o = 0.5 GeV, 
D = \ GeV'^ 04811 . In the following, we consider three quark 
flavors q = u,d,s, which are independent of each other in our 
model. We take the current-auark masses m„ j — for sim- 
plicity and Ms = 115 MeV ll48ll . The dependence of the T = 
EOS and the structure of NSs on a is discussed in Ref. [19]. 
For a w 0, the EOS of QM is too stiff to obtain a hadron-quark 
phase transition in NSs, while for aKi°°, the model is reduced 
to a free quark system at finite chemical potential. Herein, we 
take two typical values a ~2 and a = 4 in our calculation. 

At finite chemical potential and temperature, the quark 
propagator assumes a general form with rotational invariance 



one can calculate the quark number density and correspond- 
ingly ob tain the pressure and EOS for QM at zero temperature 
0191143114611 . In the case T > 0, we express the single-flavor 
quark number density as 

P,(M„r) = gjA^f^(\p\;^^J), (20) 

f,{\p\;^^,T) - y £ ti-Dh7A(p,«„;Ai^,r)],(21) 

where g = 2Nc = 6 is the quark degeneracy and the trace is 
over spinor indices only. Then, due to asymptotic freedom at 
large chemical potential /Xuv, the quark entropy density and 
the difference of pressure at jUuv can be approximated as 
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where s^^^^ {nuy ,T) is the quark entropy density in a free 
quark system. Therefore, if asymptotic freedom is approached 
continuously, which is realized in the Wigner phase in our 
model with a > 0.3, the single-flavor quark thermodynamic 
pressure at finite chemical potential and temperature can be 
obtained as 
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The total pressure for the quark phase is given by summing 
contributions from all flavors. For comparison with the bag 
model, we write the pressure of the quark phase as 

PQ{lXu,li,i,lis,T)= ^ Pg{ii,J)-Bus, (25) 

q=u.d,s 

with the bag constant 

Bds = - E PM.0,0). (26) 

and the reduced pressure Pg{Hy,T) = Pgin,,,!) -Pq{Hqfi,Q) 
from the integrals of quark number density and entropy den- 
sity. As in Ref. IJ^, we set Bds = 90 MeV with ^„.o = ji^.o = 
and /Xj as the value of the starting point of the deconfined 
phase of strange quarks. We choose /luv = 1 GeV in our cal- 
culation, which is high enough as can be seen in the upper 
panels of Figs.[T]and|2] 

Correspondingly, the total baryon number density, entropy 
density, internal energy density, and free energy density for 
the quark phase are obtained as 

pB = 3 E P?(M<?,r), 



(27) 



Sg{p,co„;ix,T) ' = ijpA{p^,co„;iXg,T)+B{p^,co„;iiy,T) 
+ iY4{(0„ + ipiq) C{p^,(o„-iig, T) . (19) 

Here we ignore the possibilities of color superconductivity 
and other structures 1,18.1 . From the quark propagator. 
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FIG. 1: Quark number density (upper panel), pressure Pq (central 
panel), and entropy density (lower panel) for massless quarks at finite 
chemical potentials and temperatures. Different curves correspond to 
different values of the parameter a = 2 (dotted blue lines), 4 (dashed 
red lines), oo (solid black lines) and different temperatures T = (thin 
lines), r = 40 MeV (thick lines). All Suj{iiJ = 0) = 0. 



temperature effects is similar in the DSM as in the free quark 
system. The quark entropy density increases with the interac- 
tion strength parameter 1 /a of the DSM. 

In summary, both at zero and nonzero temperature, the 
quark number density and pressure Pq in the DSM are reduced 
compared to the free quark system. Since the typical temper- 
atures are much smaller than the chemical potential, the tem- 
perature effects are usually small. 
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In Figs. [T] (for massless up and down quarks) and |2] (for 
strange quarks) we show the quark number density (upper 
panel), pressure P^ (central panel), and entropy density (lower 
panel) at finite chemical potentials and temperatures. Both 
at zero and nonzero temperature, when the chemical poten- 
tial is about a few hundred MeV, the quark number density 
in our model (DS2 and DS4) is lower than that of the free 
quark system DSoo. This leads to a reduction of the pressure 
Pq. Qualitatively, the increase of quark number density due to ^ '^^y'^j 



FIG. 2: Same as Fig.[T] but for strange quarks with m,s = 1 15 MeV. 



B. The MIT bag model at finite temperature 

We review briefly the description of the bulk properties of 
uniform QM at finite temperature within the MIT bag model 
[56.1 . In its simplest form, the quarks are considered to be free 
inside a bag and the thermodynamic properties are derived 
from the Fermi distribution, i.e.. 



fq{k-jXqJ) = f^{k)-f-(k), 

fHk) = 



l+exp[{Eqik)Tfici)/T] 



(31) 
(32) 



' + k^ and jlq being the auxiliary quark 

chemical potentials. The quark number density is obtained 
from Eq. (|20| | and the energy density and free energy density 
for the quark phase are given as 
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where B is the bag constant and s^^^ the entropy density of 
a noninteracting quark gas. The 'true' quark chemical poten- 
tial, entropy density, and pressure are given by the general 
relations 
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The corresponding expressions at T = can be obtained by 
eliminating the antiparticles and substituting the particle dis- 
tribution functions by the usual step functions. We have used 
massless u and d quarks, and m^ = 150 MeV. 

It has been found |57] that within the MIT bag model (with- 
out color superconductivity) with a constant bag parameter 
B ~ 90 MeV fm^^, the maximum mass of a NS cannot ex- 
ceed a value of about 1.6 solar masses. Indeed, the maxi- 
mum mass increases as the value of B decreases, but too small 
values of B are incompatible with a hadron-quark transition 
at baryon density pg > 2-3 po in nearly symmetric nuclear 
matter, as demanded by heavy-ion collision phenomenology. 
(These baryon densities are usually reached in numerical sim- 
ulations [58] of heavy ion collisions at intermediate energies 
without yielding indications of 'exotic' physics.) 

In order to overcome these problems, we have introduced 
in a phenomenological manner a density-dependent bag pa- 
rameter B{pb) [57]. This allows one to lower the value of 
B at large density, providing a stiffer QM EOS and increas- 
ing the value of the NS maximum mass, while at the same 
time still fulfilling the condition of no phase transition below 
pB ~ 3po in symmetric matter In the following we present 
results based on the MIT model using both a constant value 
of the bag parameter, and a Gaussian parametrization for the 
density dependence. 
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withBoc = 50 MeV fm^^ Bq = 400 MeV fm-^ and j3 = 0.17, 
see Ref. IdTII . In this paper, we disregard possible dependen- 
cies of the bag parameter on temperature and individual quark 
densities. For a more extensive discussion of this topic, the 
reader is referred to Refs. (|57| |59|. 



IV. RESULTS AND DISCUSSION 
A. Matter in beta equilibrium 

In neutrino-trapped beta-stable nuclear matter the chemical 
potential of any particle i — n,p,A,l,^,u,d,s,l is uniquely de- 
termined by the conserved quantities baryon number B,, elec- 
tric charge Qi, and weak charges (lepton numbers) L^ , L^ 
with corresponding chemical potentials plbt^Iq, jJ-ie , I-LLu • 



where B,-, Qi,Lf ,L: are the corresponding charges of each 
particle. The relations of chemical potentials and densities for 
hadrons and quarks are given in the previous sections, and we 
treat leptons as free fermions. With such relations, the whole 
system in each phase can be solved for a given baryon density 
Pb = HiBrPi, imposing the charge neutrahty condition 



Eap'-o 



and lepton number conservation 
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We fix the lepton fractions to Y^ = 0.4 for neutrino-trapped 
matter at T = 40 MeV, and we neglect muons and muon- 
neutrinos due to their low fractions, hence Y^ =Q . When 
the neutrinos Ve are untrapped, the lepton number is not con- 
served any more, the density and the chemical potential of Ve 
vanish, and the above equations simplify accordingly. 

In Fig. [3] we display the relations between baryon chemi- 
cal potential, baryon density, and pressure for the pure hadron 
and quark phases of beta-stable matter Thin curves repre- 
sent the hadronic phases obtained with the Argonne Vi% NN 
potential supplemented by microscopic TBF (solid black) or 
phenomenological Urbana TBF (dotted black), and including 
also hyperons (dashed black). The thick curves are the re- 
sults for beta-stable QM, i.e., the Dyson-Schwinger results 
with a — 2 (dash-dot-dotted blue) and a = 4 (dotted red), 
whereas the other two lines represent the results obtained with 
the MIT bag model, using either a constant bag parameter 
B = 90 MeV fm^^ (dashed green) or a density-dependent 
B(p) (dash-dotted pink). The left-hand panels display results 
for cold beta-stable matter, the central panels for hot (T = 40 
MeV) and untrapped matter, and the right-hand panels for hot 
and trapped matter The crossing points of the baryon and 
quark pressure curves represent the transition between baryon 
and QM phases under the Maxwell construction. The projec- 
tions of these points (vertical lines) on the baryon and quark 
density curves in the upper panels indicate the corresponding 
transition densities from low-density baryonic matter to high- 
density QM. 

We notice that the phase transition from hadronic to QM oc- 
curs at low values of the baryon chemical potential when the 
MIT bag model is used to describe the quark phase. This holds 
true irrespective of the EOS adopted for the hadronic matter, 
since the TBF effects are only important at high densities. On 
the contrary, with the DSM for QM the phase transition oc- 
curs at higher values of the density, since the EOS is generally 
stiffer than the hadronic one. In this case, the onset and the 
width of the density range characterizing the phase transition 
depend strongly on the EOS used for the hadronic phase. 

Thermal effects play an important role, since they shift the 
onset of the phase transition to lower values of the baryon 
density. Neutrino trapping lowers even more the transition 
density in the MIT case, but increases it for the DSM. This 
is because of the different behavior of the nucleonic EOS at 
low and high densities: At low density, corresponding to the 
transition density in the MIT case, neutrino trapping increases 
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FIG. 3: (Color online) Baryon density (upper panels) and pressure (lower panels) as function of the baryon chemical potential for cold 
untrapped matter (left panels), untrapped matter at T = 40 MeV (central panels), and trapped matter at T = 40 MeV (right panels). See text 
for details. 



the pressure and decreases the transition density. At large den- 
sity, corresponding to the transition density in the DSM case, 
nuclear matter has a larger symmetry energy, and the decrease 
of nucleon pressure overwhelms the increase from the leptons. 
Therefore, neutrino trapping decreases the total pressure and 
increases the transition density. We also notice that with the 
DSM no phase transition exists if the hadronic phase contains 
hyperons and the EOS is very soft, or if the parameter a is 
too small and the EOS of QM is very stiff, in analogy with the 
zero-temperature case IITqII . 



B. The hadron-quark phase transition 

We found in Ref. [19] that hybrid NS with the DSM ob- 
tained using the Maxwell construction are unstable. We there- 
fore study in the following the phase transition at finite tem- 
perature with the Gibbs construction [20., 60, 611 , which deter- 
mines a range of baryon densities where both phases coexist, 
yielding an EOS containing a pure hadronic phase, a mixed 
phase, and a pure QM region. In the Gibbs construction, as 
suggested by Glendenning ||60|] , the stellar matter is treated 
as a two-component system, and therefore is parametrized by 
two chemical potentials. Usually the pair (\ig.\i^, i.e., elec- 
tron and baryon chemical potential, is chosen. By imposing 
mechanical, chemical, and thermal equilibrium, it turns out 



that the pressure is a monotonically increasing function of the 
density, at variance with the Maxwell construction, where a 
plateau in the pressure vs. density plane exists, thus character- 
izing the phase transition. We note that the Gibbs treatment 
is the zero surface tension limit of the calculations including 
finite-size effects [61, 62]. The Gibbs phase transition has 
been widely studied in the literature, and the formalism will 
not be repeated here. 

In Fig.|4]we display results for the EOS involving the Gibbs 
hadron-quark phase transition, comparing calculations us- 
ing phenomenological (upper panel) and microscopic (lower 
panel) TBF. We adopt similar conventions as in Fig.|3]for the 
QM EOS, but with thin lines indicating the pure hadron or 
quark phases and thick lines the mixed phase region. Similar 
as in Fig. [3] for the phase transition with the Maxwell con- 
struction, thermal effects alone (central panels) shift the onset 
of the mixed phase to smaller density, while neutrino trap- 
ping delays the onset. Taking all effects into account, the EOS 
comprising the phase transition turns out to be softer than the 
hadronic one. 

Comparing the phase transition constructed with the DSM 
and that obtained using the MIT bag model, they turn out to be 
quite different, as predicted from Fig. |3] In the former case, 
the onset of the phase transition is shifted to a larger bary- 
onic density and the mixed phase is extended to a much larger 
region. Different nuclear EOSs also affect the onset and the 
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FIG. 4: (Color online) Pressure vs. baryon density of NS matter with the Gibbs phase transition construction for different models. Phenomeno- 
logical (upper panels) or microscopic (lower panels) TBF are used for the hadronic EOS. The pure hadron or quark phases are marked with 
thin curves and the mixed phase region with thick curves. The pure quark phase regions with the DSM in the upper panel and with DS2 in the 
lower panel lie beyond the range we plot. 



width of the mixed phase, in particular when the DS model is 
used. In the case of a stiff EOS for the hadronic phase with 
the microscopic TBF, the onset of the mixed phase is at lower 
density and the mixed phase has a smaller region, compared 
to the case of phenomenological TBF. These all have conse- 
quences for the internal structure of a PNS, as discussed in the 
next subsection. 



C. Proto-neutron star structure 

The stable configurations of a (P)NS can be obtained from 
the well-known hydrostatic equilibrium equations of Tolman, 
Oppenheimer, and Volkov yj for pressure p{r) and enclosed 

mass m{r), 



dp 
dr 
dm 
dr 



Gme (l +ple){l +47:r^p/m) 



1 — 2Gm/r 



= Anre, 



(42) 
(43) 



once the EOS p{e) is specified (G is the gravitational con- 
stant). For a given central value of the energy density, the nu- 
merical integration of Eqs. (l42t and (|43] i provides the mass- 
central density relation. In the low-density range, where 



nucleonic clustering sets in, we cannot use the BHF ap- 
proach, and therefore we join lUlll the BHF EOS to the finite- 
temperature Shen EOS 06311 . which is more appropriate at den- 
sities below p < 0.07 fm^ , since it does include the treatment 
of finite nuclei. 

Our results, using the different EOSs introduced in the pre- 
vious section, are displayed in Fig. |5] which shows the grav- 
itational mass (in units of the solar mass M© = 1.98 x lO^^g) 
as a function of the central baryon density. We use the same 
conventions as in the previous figures: The black solid curves 
represent the calculations performed for purely hadronic mat- 
ter, using either phenomenological (upper panels) or micro- 
scopic (lower panels) TBF. The colored broken curves denote 
stellar configurations of hybrid (P)NSs with neutrino-free (left 
and central panels) and neutrino-trapped matter (right panels). 

We notice that the value of the maximum mass decreases in 
neutrino-free matter due to thermal effects, both in the purely 
hadronic case and including the hadron-quark phase transi- 
tion, where the decrease depends on the EOS used for the 
quark phase and turns out to be more relevant for the DS 
model. On the contrary, neutrino trapping decreases further 
the maximum mass of purely hadronic stars, but increases the 
value of the maximum mass of hybrid stars, overcoming the 
thermal effect. Therefore in this case a delayed collapse sce- 
nario is possible, just as for hyperon stars [,12,1. 
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FIG. 5: (Color online) NS gravitational mass vs. central baryon density. Phenomenological (upper panels) or microscopic (lower panels) TBF 
are used for the hadronic EOS. 



One also notices a dependence on the EOS used for the 
hadronic phase, which is more important for the DS model, 
where the QM onset takes place in a range of densities where 
TBF play an important stiffening role, and this explains the 
different values of the maximum mass. On the contrary, with 
the MIT bag model the transition takes place at low baryon 
densities where the different TBF behave similarly, hence the 
small influence. The proper values of the maximum mass 
are summarized in Table U and a maximum mass exceeding 
2 Mq for the hybrid configurations can be obtained with the 
microscopic TBF and suitable values of the parameter a in 
the DSM. 



V. CONCLUSIONS 



TABLE I: Properties of (P)NS maximum mass configurations. Mi- 
croscopic TBF are used. 



Extending the work of Ref. 111911 to finite temperature, 
we have studied the properties of hybrid PNSs based on a 
hadronic BHF EOS involving different nuclear TBF joined 
via a Gibbs construction to the DSM quark EOS. This is only 
possible if hyperons are excluded from the hadronic phase. 

We find a sizeable dependence of the results, in particu- 
lar for the PNS maximum mass, on the nuclear TBF and on 
the quark model employed. The DSM model allows to reach 
larger masses than the MIT model, depending on the value of 
the interaction strength parameter. 

In all cases finite temperature reduces slightly the maxi- 







M/Mq 


R(km) 


Pc/PO 




N (mTBF) 


2.30 


10.85 


5.86 


untrapped 


DS2 


1.87 


11.02 


6.54 


T=0 


DS4 


1.60 


10.61 


7.09 




MIT, B=90 


1.51 


9.02 


9.27 




MIT B=B(p) 


1.53 


8.37 


10.23 




N (mTBF) 


2.28 


13.14 


5.72 


untrapped 


DS2 


1.78 


14.68 


6.47 


T=40 MeV 


DS4 


1.53 


15.46 


6.65 




MIT B=90 


1.51 


11.62 


8.87 




MIT, B=B(p) 


1.51 


10.37 


10.13 




N (mTBF) 


2.19 


12.85 


6.18 


trapped 


DS2 


2.05 


15.33 


5.70 


T=40 MeV 


DS4 


1.89 


17.55 


5.26 




MIT, B=90 


1.61 


15.82 


6.98 




MIT, B=B(p) 


1.55 


13.37 


8.76 



mum mass, but neutrino trapping increases it sufficiently for 
hybrid stars, such that a delayed coUapse might be possible in 
principle. 
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